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Abstract. The paper is concerned with the Barenblatt-Biott model in the 
theory of poroelasticity. We derive a guaranteed estimate of the difference 
between exact and approximate solutions expressed in a combined norm that 
encompasses errors for the pressure fields computed from the diffusion part 
of the model and errors related to stresses (strains) of the elastic part. Es- 
timates do not contain generic (mesh-dependent) constants and are valid for 
any conforming approximation of pressure and stress fields. 



1. Introduction 

The standard mathematical model for diffusive flow in an elastic porous media 
is the Biot's diffusion-deformation model of poroelasticity [3] based on coupling 
between the pore-fluid potential and the solid stress fields. The basic constitutive 
equations relate the total stress to both the effective stress given by the strain of 
the structure and to the potential arising from the pore-fluid. The model consists 
of a momentum balance equation combined with Hooke's law for elastic deforma- 
tion, and a continuity equations combined with Darcy's law. Originally, Biot's 
model was designed for homogeneous porous media or single porosity media. The 
representation of porosity and permeability in naturally occurring materials often 
requires several distinct spatial scales. As for instance, in reservoir model, the pres- 
ence of heterogeneities like highly permeable channels has a significant impact on 
the flow properties of reservoir rock. Two of more scales of permeability are usually 
observed, which is also referred to as dual permeability models. 

Studies suggest that even for single-phase flow in relatively simple porous media, 
such as sandstone, the fluid flows through a very small portion of the pore space, 
while a greater part of it remains stagnant. A connected system of highly permeable 
channels, characterized by relatively simple pore space geometry, provides fluid flow 
through the reservoir. The remainder of the reservoir, characterized by tortuous 
pores and pore throats, is significantly less permeable. The highly permeable chan- 
nel component of a reservoir is relatively small, and the remainder of the reservoir 
contains most of the fluid. This contrast leads to the dual medium model of reser- 
voir rock, originally proposed by Barenblatt et al. [I] in the rigid case. According 
to this model, the fluid flow in matrix blocks is local, and only the local exchange 
of fluid between individual blocks and the surrounding high permeable channels is 
supported. This model contains a system of two diffusion equations, one for each 
component, coupled by a distributed exchange term that, in its simplest form, is 
proportional to the difference in potential between fluids in the two components. 
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A combination of the Barenblatt's double-diffusion approach and Biot's diffusion- 
deformation theory leads to what we call the Barenblatt-Biot poroelastic model 
representing double diffusion in elastic porous media. It takes the form 

- V • (Le(u)) +aiVpi +a 2 Vp 2 = f(x,t), 

(1.1) cypx ~ V • (feiVpi) +«iV • ii + «(pi -p%) = hi(x,t), 
C2P2 - V • (fc 2 Vp 2 ) + a 2 V • u + k(p 2 -Pi) = h 2 (x,t), 

u is the displacement of the solid skeleton and p\ and p 2 are the fluid potentials in 
the respective components. With the vector gradient operator V, the linear Green 
strain tensor e(-) writes 

(1.2) £ ( U ):=i(Vu+(Vuf). 

The fourth-order elastic stiffness tensor L defines a stress tensor er using the Hook's 
law 

a := Le(u). 

In general, the permeabilities k\ and fc 2 may be heterogeneous and anisotropic 
tensors, which may be functions of the deformation. Herein, we will neglect this 
dependence and only consider constant, scalar and homogeneous permeabilities. 
Constants ot\ and a 2 measure changes of porosities due to an applied volumetric 
strain. Mathematical analysis of this model based on the theory of implicit evolution 
equations in Hilbert spaces is elaborated in [T2] , 

We note that multiple continua models are applicable to several other porous 
media problems. We mention two cases in particular. Firstly, contaminant trans- 
port experiments clearly indicate that the particle dispersion is non-Fickian, as 
reviewed in [2]. This makes both dual and multiple continua models of interest, 
with dual media approaches already common in applications. The use of more than 
two flowing continua was argued by Gwo et al. [5], while a single flowing continua 
coupled to multiple non- flowing continua (traps) is also reviewed in [2] . The second 
application is heat transfer in fractured rocks, in particular related to modelling of 
geothermal heat extraction. Here, the slow interaction of diffusive heat transfer in 
rock has to be modelled together with fast fluid flow in fractures. The state of the 
art approach is to use multiple continua, frequently as many as four or more [7] . 

Our focus in this paper is to derive guaranteed and computable bounds of ap- 
proximation errors the static Barenblatt-Biot system 

- V- (Le(u)) + a 1 Vp 1 + a 2 Vp 2 = f(x), 

(1.3) -V-(*iVpi) + /s(pi-pa) = hi(x), 
-V ■ (fc 2 Vp 2 ) + k(p2 -pi) = h 2 (x), 

which is considered in bounded connected domain ft C lZ d with Lipschitz continu- 
ous boundary T. There are various boundary conditions motivated by hydrological 
applications, among which four boundary conditions, applicable to different parts 
of the boundary r = (J 1^ represent the most typical cases. 
1. Saturated land surface, T\, with infiltration and evaporation is modelled as 

(1.4) er(u)n = (normal stress free condition), 

(1.5) VTi = n • (— fciVpi — /c 2 Vp 2 ) (normal fluid flux), 
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where n is the unit outward normal vector and ipr 1 is a given function. We complete 
this boundary condition by specifying that the normal component of the potential 
gradients at the boundary are equal 

(1.6) n-V(pi-p2)=0. 

In the case of constant k (that we consider in this paper), the condition (| 1 . 5[) reads 

(1.7) VTi = -( fc i + fc 2)n • Vpi = -(h + fc 2 )n • Vp 2 , 
which is in fact a version of the Darcy law at the boundary. 

2. Boundary to sea with a constant fluid potential (we call this boundary V^)- Here 
we also impose normal stress as in (|1.4[) . However, the boundary conditions for the 
potentials are of the Dirichlet type, i.e., 

(1-8) p 1 =p 2 =pr 2 - 

3. Internal boundary with known head (T 3/ ). This may represent either a fixed 
potential pumping well or the potential at some measurement point. We model 
this as a no displacement boundary with Dirichlet conditions for the potentials as 
at T 3 , i.e., 

(1.9) u = 0, 

(1.10) p 1 =p 2 =pr :i - 

4- Impermeable bedrock, T 4 . Here, we impose no displacement (as for T 3 ), and zero 
normal flux (as for Equation (|1.5I) with i/jr 4 = 0). 

For the unique solvability of the diffusion problem one has to assume that 

meas(r 2 U T 3 ) ^ 0. 

2. VARIATIONAL FORMULATION OF THE DOUBLE DIFFUSION SYSTEM 

Since the displacement u is only involved in the first equation of system (|1.3[) , a 
double-diffusion problem 

(2.1) - V- (AiVpi) + «(pi -P2) =hi(x), 

(2.2) -V ■{k 2 Vp2) + n{p2-pi) = h 2 {x) 

is studied separately. It describes the steady flow of slightly compressible fluid in 
a general heterogeneous medium consisting of two components. Henceforth, we 
consider this problem with the Dirichlet boundary conditions pi — P2 — Pv on T. 
Let p be a function with square summable coefficients that satisfies this boundary 
condition. It is convenient to rewrite the problem in terms of new functions 

Pi -=Pi -P, P2 -=P2 - p. 
Then, a weak formulation of (|2.1[) - (|2.2[) leads to 

Problem 1. Assume that (hi,h 2 ) € L 2 (f7,^ 2 ). Find p = (pi,p 2 ) € H^n^U 2 ), 
satisfying the system of variational equalities 

J feiVpi • Vqi + J ft(pi - p 2 )qi dx = J (hi(x)qi - fciVp • Vqi) dx 

an n 

J k 2 Vp 2 ■ Vq 2 + J ft(p 2 - pi)q 2 dx = J (h 2 {x)q 2 - k 2 Vp ■ Vq 2 ) dx 



(2.3) 



Q 



for all testing functions q— (qi,q 2 ) G Hq(£1,TZ 



2 > 
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This problem can be represented in a general form (which also encompasses 
other, more complicated models of porous media). For this purpose, we introduce 
the spaces 

(2.4) Q := H*(Cl,H*), Y := L 2 (fi, lZ 2d ), 
and the corresponding dual spaces 

(2.5) Q*:=iT- 1 (n,ft 3 ), Y* := L 2 {n,K 2d ). 

Hereafter L 2 norms of all functions in Q are denoted by ||-||q. Duality pairings of 
(Q, Q*) and (Y, Y*) are denoted by (•, •) and ((•, •)), respectively. Also, we introduce 
a bounded linear operator A e C(Q,Y) and its adjoint operator A* e C(Y*,Q*) 
by the relations 

(2.6) Aq:= (V qi ,Vq 2 ), A*Y* = (- divj/J, - divy*) T . 

The operators A and A* satisfy the relation representing integration by parts 

((¥*, Aq)) = (A*Y*,q) for all Y* £ F*,q£ Q, 
which can be written componentwise as 

(2.7) J (Y*-Vqi+Y*-Vq2)da; = -y (qi div Y* + q 2 div Y*) dx, 

n n 

where q = (qi,q 2 ) and Y* = (Y*, Y|). Now Problem Q] can be represented in the 
form: Find p 6 Q such that the equality 

(2.8) a(p,q) = J(<l) 

holds for all q S Q. The bilinear form <z(-, •) and the linear form /(•) arc dchned as 

a(p, q) := J (Ap : (AAq) + p • Bq) dx, 

n 

i(q) := y (h • q - CAq) dx, 

n 

A, B and C are matrices formed by material dependent constants fci, &2, &3, 

ki 0\ / k -k\ r ._ /feiVp 

fc 2 y ' y — k k J ^ — { yfc 2 Vp 

and h is the right hand side vector 

M 

Remark 1. We note that the symmetric matrix A is a positive definite matrix iff 
fci and fc 2 are positive (since A£ • £ > min{fci,fc 2 } ||£|| 2 for all £ 6 However, 
B is symmetric but only positive semi-definite in case of the positive parameter k, 
and its one-dimensional kernel is generated by the vector (1, 1) T . 

Remark 2. If p is sufficiently regular (so that A*C belongs to Y*), then 

2(q) := y"( h -q- A*Cq)d.T = J h ■ qdx, 
n n 
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where 



hi — div fci Vp 
hi — div &2 



It is easy to verify that (|2.8j) is the necessary condition for the minimizer of the 
following convex variational problem. 

Problem 2. Find p G Q satisfying 

(2.9) F(p) + G(Ap) = inf {F(q) + G(Aq)}, 

w/iere 

(2.10) F:Q~>K 7 F(q) := i / q • Bqdx - Z(q), 



(2.11) G : F -> ft, G(Aq) :=]^ J kq : (AAq) Ac. 

Theorem 1 (existence of unique solution). Assume that ki,k2 > and re > 0. 
Then, there exists a unique solution p £ Q o/ Problem which also represents the 
solution of ProblemUi 

Proof. Existence of the unique minimizer follows from known results in the calculus 
of variations. Indeed, under the give assumptions, the functional F(-) + G(A-) is 
strictly convex and coercive in the reflexive space Q. □ 

3. A POSTERIORI ERROR ESTIMATE OF THE DOUBLE DIFFUSION SYSTEM 

In this section, we derive guaranteed and directly computable bounds of the 
difference between exact and approximate solutions. Our analysis is based upon a 
posteriori error estimation methods suggested in [6l [10] . Following the chapters 6 
and 7 in [6], first we need to find explicit forms of dual functionals 

F* : Q* F*(A*Y*) := sup{(A*Y*, q) - F(q)}, 

(3 I) qGQ 

[ ' ' G* : y* — > 1Z, G*(Y*) := sup {((Y*, Aq)) - G(Aq)}. 

Aqer 

and the corresponding compound functionals 

D F : Q x Q* D F (q,A*Y*) := F(q) + F*(A*Y*) - (A*Y*,q) , 

^ 3 ' 2 ^ D G :YxY*^TZ, D G (Aq, Y*) := G(Aq) + G*(Y*) — ((Y*, Aq)) . 

By the the sum of Dp and Dq, we obtain the functional error majorant 

(3.3) M(q, Y*) := D F (q, A*Y*) + D G (Aq, Y*), 
which provides a guaranteed upper bound of the error: 

(3.4) ia(p-q,p-q) < Af(q,Y*) for all Y* E Y* . 

The majorant is fully computable and depends only on the approximation qGQ 
and arbitrary variable Y* G Y*. 
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Lemma 1 (dual functionals) . For fci, > and n > 0, it holds 

(3.5) G*(Y*) = ^ J A _1 Y* : Y* dx, 

f J (A*Y* + h) 2 dx ifA*y* 1 +h 1 +A*yZ + h 2 =0, 

(3.6) F* (A* Y*) = \ n 

[ +oo otherwise. 

Proof. The derivation of G* (Y* ) is straightforward, see [6] . The singularity of the 
matrix B makes the computation of F*(A*Y*) more technical. 

F*(A*Y*) = sup{(q,A*Y*)-F(q)} 
qeQ 

> sup {(q,A*Y*)-F(q)} 

q£Q:qi=q 2 

= sup {{ qi ,A*Yt+A*Y* 2 )-F( qi , qi )} 
qie-Ffi(n) 

sup {(q 1 ,A*Y* 1 +h 1 +A*Y*+h 2 )} 

if A*Yl +hi+ A*Y* 2 + h 2 = 0, 
+oo otherwise. 

Thus, finite values _F*(A*Y*) are attained only on the subspace 

(3.7) A*Y* + hi + A*Y; + h 2 =0 
and we must specially consider this case. It holds 

F*(A*Y*) = sup{(q,A*Y*) -F(q)} = sup{(q,A*Y* +h) - J /Bq-qdx} 
qeQ qeQ 2 J 

n 

(use the constrain A*Y^ + h 2 = -(A*Y* + hi)) 
= sup {(qi - q 2 , A*Y* + hi) - - / «(qi - q 2 ) 2 dx} 

(qi,q 2 )eQ Z J 

(supremum is obtained for qi — q 2 = (A*Y* + hi)/n) 

= i- J {A*Yl + hi) 2 dx = ^J [(A*Y* + h^ 2 + (A*Y 2 + h 2 ) 2 } dx 
n Q 

= — [ (A*Y* + h) 2 dx. 
4k J 



□ 



Remark 3. We note that (|3.7[) is a weaker restriction than the sum of two equilib- 
rium relations A* Y^ + hi = and A* Y| + hi = 0, which one would await from the 
general theory. In other words, our analysis shows that the strict equilibrium of 
the dual variables in the componentwise sense is not required in the couple system. 

After the substitution of p. 51) and Q3.6P in the definition (|3.2I) , we obtain explicit 
expressions for the compound functionals. 
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Lemma 2 (compound functionals) . It holds 

(3.8) L> G (Aq,Y*) = iy*A(Aq- A -1 Y*) : (Aq - A _1 Y*)da;, 



| / Bq • q dx + i /(A* Y* + h) 2 da; 

(3.9) L» F (q,A*Y*) = <j n i/ A*Y* + hi + A* Y?j + h 2 = 0, 

foo otherwise. 

According to (|3.23|) . the sharpest bound of a(p q, p — q) is provided by the 
estimate 

(3.10) ia(p-q,p-q)< y inf B M(q,Y*). 

Since M(q, Y*) = +oo if Y* does not satisfy (|3.7[) . we must restrict ourselves to 
arguments Y* £ 1^, where 

(3.11) 1£ := G Y* : A*YJ + + A*Y* 2 + h 2 = a.e. in 0}. 

To construct an element of Yfi requires an exact equilibration procedure, which 
have been studied for a Poisson problem in [3] . Below, we show a way to avoid the 
constrain (|3.11[) by a special penalty term added to the functional major ant. We 
define 

(3.12) Y£ v := {(Y*,Y*) e Y* : A* Y* + A* Y* e L 2 {fl)} 

and note that Y h * C Y£ iv (since hi,h 2 £ L 2 (Q)). Further we decompose 

Y* = Y* + (Y* — Y*) 

with Y* £ Yf iv and we extend the dual functionals Dq and Dp by the new variable 
Y*. We rewrite (O as 



L> G (Aq, Y*) = ~ J A(Aq - A^Y*) : (Aq - A^Y*) dx + 
a 

+ J (Aq- A _1 Y*) : (Y* -Y*)dx+^ J A _1 (Y* - Y*) : (Y* -Y*)dx 
n a 

and use the inequality 2Mi : M 2 < $\M\ : M x + ^-M 2 : M 2 valid for all matrices 
Mi , M2 and for all /3% > to bound the middle term as 

(3.13) (Aq - A _1 Y*) : (Y* - Y*) = A 1/2 (Aq - A _1 Y*) : A~ 1/2 (Y* - Y*) 

< — A(Aq - A" 1 **) : (Aq - A _1 Y*) + -^-A^Y* - Y*) : (Y* - Y*). 
2 2pi 

Obviously, the middle terms adds to the left and the right terms in Z) G (Aq, Y*) 
above and the modified compound functional reads 

L> G (Aq,Y*,Y*) :=^r^- J A(Aq - A _1 Y*) : (Aq - A _1 Y*) dx 

^ 1 ° 1 f 

+ ( 2 + 2p\ } J A_1 ( Y *- Y *) : ( Y *- Y *)^- 



s 
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It also contains a scalar factor j3\ > that value can be chosen arbitrarily. Similar 
technique is used to modify the compound functional D F {q, A*Y*). For the second 
integral in (|3.9I) . we have 

J (A*Y* + h) 2 dx < (1 + I3 2 ) J (A*Y* + h) 2 dx + (1 + 1.) J (A*(Y* - Y*)) 2 dx, 
n n n 

where /?2 > 0. Therefore, a modified dual functional reads 

D i r(q,A*Y*,A*Y*) :=1 / Bq • qdz + ^-(1 + ft) / (A*Y* + h) 2 dx 
n n 

+ 4^ (1 + ^ } /( A *( Y *-Y*)) 2 ^. 



(3.15) 



S2 



By adding (|3.14l) and (j3. 15[) . we extend the functional majorant (13.31) to 

(3.16) M(q,Y*,Y*) := D F (q, A*Y*, A*Y*) + L> G (Aq, Y*, Y*), 
in which arbitrary variables satisfy the constrain 

(Y*,Y*)ey,;xy;„,. 

Clearly, the original and extended majorants satisfy the inequality 

(3.17) io(p - q, p - q) < M(q, Y*) < M(q, Y*, Y*) 

for all Y* G > 0,/?2 > 0. This estimate is sharp in the sense that there are 

no irremovable gaps in the inequalities. Indeed, if we set Y* = Y* = Ap and tend 
/?i and fa to zero, then M(q, Y*,Y*) tends to M(q,Y*) (and even to the exact 
error |a(p — q, p — q), cf. (|3.10p ). 

3.1. An upper estimate of M(q,Y*,Y*). Let us denote Y* = (YJ,Y|) and 
Y* = (Y* , Y2 ) and consider a particular subspace 

(3.18) (Y*, Y*) £ {Y£ x Y d * m : A*Y? + h t = 0, Y 2 * = Y 2 * a.e. in Ci}. 
In this subspace, it holds (cf. §'2.6]i ) 

J(A*(Y* - Y*)) 2 dx = y*(div(Y* - YDfdx = J {divYl - In) 2 dx. 
an n 
Therefore, D F (q, A*Y* , A*Y*) defined in (pT5]) simplifies as Y* -independent 

(3.19) D F (q,A*Y*) := i / Bq • qdx + i-(l + /3 2 ) / (A*Y* + h) 2 dx 

n a 

+^(l+^)/(divY 1 *-/ ll ) 2 ^ 
o 

and only Y*-dependent functional in D G (Aq, Y*, Y*) defined in (|3.14p writes 

(3.20) / A _1 (Y* - Y*) : (Y* - Y*) dx = J fcf^Y* - Y*) ■ (Y* - Y*) dx. 
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Lemma 3. Let us define a space 



Y hl := {Y* G L 2 (Cl) d : A*Y* + h x = a.e. in CI}. 



Then, for all Y| G iJ(div; CI), it holds 



inf 



Y* -Yl 



dx < C z 



divY* + /n 



where C > satisfies Friedrichs' inequality \\w\\ L2 ^ < C||^ w llL 2 (n) valid for all 

Proof. It follow from Theorem 6.1 from [TT] by the modification related to the fact 
the we consider vector arguments. □ 

Application of Lemma [3] to (|3.20p and the back substitution to p. 141) defines a 
Y*-independent dual functional 



(3.21) £ G (Aq,Y*):= 



1 + 01 



A(Aq - A -1 Y*) : (Aq - A _1 Y*) dx 



1 1 



+ Kl( 2 + 2fi )C2 d[v ^ +hl 
which provides an upper estimate of the quantity 

inf D G (Aq,Y*,Y*). 

Therefore, the sum of (|3.19l) and (|3.21[) defines a Y*-independent functional 

(3.22) M fti/32 (q,Y*) :=D F (q,A*Y*)+D G (A< 1 ,Y*) 

that serves as an upper bound of M(q, Y*, Y*) and provides a computable estimate 

(3.23) ~a(p - q,p - q) < Mp lt p 2 (q, Y*) for all Y* G Y* iv . 

Remark 4 (symmetric form of Dq). If we replace the subspace (|3.18p by 

(3.24) (Y*,Y*) G {Y£ x Yj iv : A*Y* + h 2 = 0,Y* = Y* a.e. in CI}, 
then, instead of (|3.21|) . we obtain 



(3.25) # G (Aq,Y*):= 



A(Aq - A -1 Y*) : (Aq - A -1 Y*) dx 



+ ^ 1 4 + ^)C 2 ||divY*+/ l2 



4. A POSTERIORI ERROR ESTIMATE FOR APPROXIMATIONS OF THE COUPLED 

SYSTEM (1.1) 

Assume that the fluid pressures pi and p 2 are resolved exactly and substituted 
to the elasticity equation (cf. 

- V • (Le(u)) = f(ar, t) + aiVpi + a 2 Vp 2 . 
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Let v be an approximation of u (this problem is considered in the same domain 
as the problem (|2.1|l - (|2.2|l ). We define the Dirichlet boundary condition by a 
function Uo G H 1 (Q,;TL d ) and assume 

v &MQ+Hl{Q,;n d ). 

Lemma 4. For every function t e Q := {a e L 2 (^;K^) : diva £ L 2 (£!;7?. d )} it 
holds 

(4.1) ||e(u - v)|| L . n < ||e(v) -^t^ + C ||divr + f - oiVpi - a 2 Vp 2 \\ n , 
where the constant C > satisfies an inequality 

(4.2) l|w|| n <C||e(w)|| L . n for all w e H^(il;TZ d ). 
and the norm \\-\\ is defined as 

\\ £ \\l-,n '■= I Le :edx - 



Proof. Estimates in chapter 6.5 in |10| which are applied to the linear elasticity 
problem with the right-hand side f — aiVpi — a2Vp2- The existence of constant C 
follows from Korn's and Friedrichs' inequalities. □ 



Remark 5. The estimate (|4.1|) is sharp with respect to parameter r. Indeed, the 
choice t = Le(it) satisfies the equilibrium condition 

(4.3) div t + f = ct\ Vpi + a-2 Vp2 

and reduces therefore (|4.1|) to the equality. 



Let qi and q2 be approximation of exact pressure fields p\ and p 2 respectively. 
By triangle inequalities, we obtain 

(4.4) 1 1 div r + f - a-iVpi - a 2 Vp 2 || < ||divr + f - aiVqi - a2Vq 2 || n 

+ l|V(pi-qi)|| n + ||V(p 2 -q2)|| n . 



Use (|4.4p and square both parts of (|4.1I) to obtain 

||e(u-v)||*. n < (||e(v)-L- 1 r|| L . n 
(4.5) +C||divr + f-aiVqi -a 2 Vq 2 || n 

+ C||V(pi-qi)|| n +C]]V(p2-q 2 )]| fi 

By the algebraic inequality 



(a + b + c) 2 < (l + /3 4 + /3 5 ) a 2 + (l + J-+/3 6 ) b 2 + (l + ^ + ^) c 2 

Pi P5 P6 
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valid for all scalars a, b, c and for all ft, ft, ft > 0, inequality (|4.5[) and the following 
inequality (ft is an arbitrary positive constant) 

(4.6) (||V( Pl -q 1 )|| + ||V(p 2 -q 2 )|| n ) 2 

< (1 + ft) ||V(pi - qi )||* + {l + h ||V(p 2 - q 2 )f n 

P3 

l + ft l + fo 
<max{ — , } a(p q,p-q 

Kl fc 2 ft 

<2max{ — Af ft:/ 3 2 (q,Y ). 

ki k 2 p 3 

Now we obtain the final estimate in terms of the coupled error norm 



(4.7) o(p - q,p - q) + || e (u - v)||£. n < (1 + ft + ft) ||e(v) - L - x t| 



L :Q 



1 



where 



1 + jr + ft ) G A ||div r + f - a x Vqi - a 2 Vq 2 ||* + 2C (q, Y*), 

P4 



C=l + C 1 + — + — max ' 



/?5 ft/ I fa feft 

This estimate holds for all r e Q, Y* e F d * w and all ft, . . . , ft > 0. 
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